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The constraints on the radion evolution in the Universal Extra Dimension (UED) model from 
Cosmic Microwave Background (CMB) and Type la supernovae (SNe la) data are studied. In the 
UED model, where both the gravity and standard model fields can propagate in the extra dimen- 
sions, the evolution of the extra dimensional volume, the radion, induces variation of fundamental 
constants. We discuss the effects of variation of the relevant constants in the context of UED for 
CMB power spectrum and SNe la data. We then use the three-year WMAP data to constrain the 
radion evolution at z ~ 1100, and the 2 a constraint on p/po (p is a function of the radion, to 



The SNe la gold sample yields a constraint on 
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be defined in the text) is [—8.8, 6.6] x 10 yr 

p/po, for redshift between and 1, to be [—4.7, 14] x 10~ ii! yr - x . Furthermore, the constraints 
from SNe la can be interpreted as bounds on the evolution QCD scale parameter, Aqcd/Aqcd,o, 



[-1.4,2.8] x 10" 



yr 



without reference to the UED model. 
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I. INTRODUCTION 



In recent years, the interest in extra dimension mod- 
els has been revived by string/M theory, in which the 
spacetime has to be multidimensional for the theory to 
be consistent [l|. If large extra dimensions exist, the 
weakness of gravity may be "explained" by the large ex- 
tra dimensional volume 0. Among the extra dimension 
models in the literature, the brane world (see 0] for a re- 
view) model seems the most popular, where the standard 
model (SM) fields are trapped in a (3+l)-dimensional hy- 
persurface, whereas gravity can freely propagate in the 
bulk. Yet another class of extra dimension models where 
all the SM fields can propagate in the extra simensions 
has also been seriously considered. In this type of model, 
so called Universal Extra Dimension (UED) model, each 
SM particle is associated with a tower of Kaluza-Klein 
(KK) particles, with the mass difference between states 
equal to the inverse of the size of the extra dimensions. 
This type of model is phenomenologically interesting be- 
cause it can be investigated in present and future collider 
experiments Q. Furthermore, the UED may explain the 
proton stability J5| , the number of generations |6], neu- 
trino masses [3, Q and fermion chirality . A cosmolog- 
ically attractive feature of UED is that its lightest KK 
partner is stable and so it is a natural candidate for dark 
matter [IE El. 



It is quite unnatural to think that the extra dimen- 
sional volume, called the radion, remains fixed per se, 
and so some mechanism must be initiated to stabilize 
the radion [l2|, [H[. However, in [3], it was shown that a 
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broad range of standard stabilization mechanisms fail to 
admit a static solution in the UED model in the matter 
dominated regime. Besides, the redshifting of matter can 
also serve as a source to drive the evolution of the radion 
[lj| ■ Thus, we may envision that the radion does not sit 
at the minimum of the potential quietly at the epochs 
of CMB formation and SNe la explosions. In the UED 
scenario, the evolution of the radion will naturally induce 
the variation of various fundamental constants. We shall 
see that, in the Einstein frame, the fermion masses, the 
gauge conplings and the QCD scale parameter Aqcd will 
acquire radion dependence. Effects of variations of extra 
dimensional volume on the fundamental constants have 
also been discussed in Ref . [H, [13] • 

By constraining the variation of constants at different 
redshifts, we can map out the allowed ranges of radion 
evolution at different epochs of the universe. Constraint 
on radion evolution in the UED scenario by Big Bang 
Nucleosynthesis (BBN) has been obtained in [l8|. In 
this paper, we seek complementary constraints from the 
Cosmic Microwave Background (CMB) and Type la su- 
pernovae (SNe la) data. 

Since the CMB anisotropics are measured quite accu- 
rately by WMAP and are sensitively dependent on 
various constants, we shall use the three-year WMAP 
data to constrain the radion evolution in UED models. 
Indeed, variations of several different fundamental con- 
stants have been constrained using CMB in the litera- 
ture. Among the fundamental constants, the variation of 
the fine structure constant a receives the most attention 
HSIHlIllIll]. Other constants constrained by CMB in- 
clude the electron mass m e [U, [2f| and the gravitational 
constant G [2(| H3, [28[. However, unification theories 
[29l l3(i| and string theory [3l| suggest that the funda- 
mental constants vary in a correlated way and should be 
constrained simultaneously. For the UED model at the 
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epoch of CMB formation, the relevant constants a, m e 
and Aqcd all depend on the radion. 

In the literature [H, [H, [H, [35| , SNe la data have been 
used to constrain the variation of the gravitational con- 
stant because the underlying physics of SNe la is dictated 
by the Chandrasekhar mass Mcb oc G~ z / 2 m~ 2 , where 
m p is the proton mass. A similar method can be applied 
to constraining radion evolution because m p depends on 
the Aqcd and hence the radion. 

However, when the radion is light enough to evolve cos- 
mologically, its mass should be comparable to the Hubble 
parameter Hq ~ 10 _34 eV. Furthermore, the coupling of 
the radion to matter is not expected to be much weaker 
than that of gravity to the matter. This means this field 
is highly constrained by tests of relativity and fifth force 
experiments [36l.l37|. In particular, the solar system tests 
have ruled out such light scalar fields in the simplest mod- 
els. A possible way out of this constraint is that the ra- 
dion may be an environment-dependent chameleon field 
[38l l39l ]. In this scenario, the effective potential of the 
scalar field depends on the density of the matter. On 
earth, the density is large and the the field is massive, 
and it can be as light as Hq in the cosmological scale. 
Because of the "thin-shell effect" of the chameleon field, 
the effective coupling is much reduced, and the solar sys- 
tem constraints can be easily satisfied. The current con- 
straints on such a light scalar field can be satisfied if the 
field is lighter than 1CT 3 eV [H. 

The rest of the paper is organized as follows. In Sec- 
tion [Til we carry out the conventional dimensional reduc- 
tion to get the radion dependence of various constants 
in the context of UED. We then present the radion de- 
pendence of various constants relevant for CMB and SNe 
la, and their effects on the CMB power spectrum and 
SNe la in Section IIIII The numerical constraints on the 
radion evolution and discussions on our constraints are 
presented in Section HVl We summarize in Section [VJ 



II. 



DIMENSIONAL REDUCTION AND LOW 
ENERGY EFFECTIVE ACTIONS 



In this section we will start from the higher dimen- 
sional action and carry out the dimensional reduction to 
derive the dependence of the constants on the radion. 
The approach here follows closely that in Ref. [HI, [4(| • 
Readers who are familiar with this type of set-up may 
go to Section IIIII directly. We shall dimensionally reduce 
the higher dimensional gravitational action to obtain its 
low energy effective action first. The higher dimensional 
action reads 



S 



1 



2k 4 , 



(1) 



where n denotes the number of compact extra dimen- 
sions, and X A represents the bulk spacetime coordinate, 
A = 0, 1, . . . , 3 + n. The higher dimensional Ricci scalar 
is denoted by ( 4+ ™)i?, G is the determinant of the full 



spacetime metric Gab , and n 2 +n is a constant related to 
the higher dimensional Planck mass. C m is the matter 
field Lagrangian density, which may include scalar fields, 
gauge fields and Dirac fermion fields. 

To proceed we shall take the metric ansatz 

ds 2 = G A BdX A dX B = g^fjcjdx^dx" + h lj (x)dy l dy j . 

(2) 

In this metric the Greek indices run over 0, 1, 2 and 3, 
while the Latin indices i and j run over the universal 
extra dimensions from 4 to 3 + n. We are only interested 
in the zero-mode of the KK expansion, and so the extra 
dimensional metric hij{x) does not depend on the extra 
dimension coordinates. Furthermore, we have assumed 
that the metric in Eq. [2] is block-diagonal because the 
vector-like connection G\, vanishes for zero-mode. The 
extra dimensions are compactified on an orbifold and the 
dimensionless coordinate y l assumes values in the interval 
[0,1]. 

After expressing ( 4 +")i? in terms of the 4D Ricci scalar 
( 4 ) R and integration by parts once, we obtain 



S = J d 4 xd n yy/=g~Vh 



+ -rg pa h ii d p h ij h kl d <r h H )+C T . 



where we have defined 



Vo 



M+n 



(3) 



(4) 



with Vo being the volume of the extra dimensions today. 

We shall work in the Einstein frame with pure Ricci 
scalar in the gravitational action. To this end we will 
apply the conformal transformation 
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with 
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(5) 



(6) 



After some algebra, we obtain the effective action 



S 



2k 2 



-~g p °Wd p h l3 h kl d a h kl 



„40 , 



(7) 



For concreteness, we further assume that the extra di- 
mensional manifold is homogeneous and isotropic. Hence 
the extra dimension metric takes the simple form: 

diag(& 2 ,6 2 ,...,6 2 ). (8) 

With the ansatz Eq. [SJ the effective action reads 



S = 



XV -9 



(9) 
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where we have defined a new scalar field, the radion a, 
as 



1 

«4 



i + 2 b n 
2n n V^' 



(10) 



We note that an effective potential has to be introduced 
to stabilize the radion. The effective potential may in- 
clude contributions from the bosonic and fermionic sec- 
tors in addition to the higher dimensional potential [13] 
or the Casimir energy [12 |. 

We now study the matter sector in more details. Let 
us begin with the minimally coupled scalar field 



'Scalar 



^G AB d A ^d B 4>-v{4>) 



(11) 

Because we are only interested in the low energy zero 
mode of the KK tower, the field cj) is independent of the 
extra-dimensional coordinates, which can be integrated 
out. After the conformal transformation Eq.O redefining 
the field as <j> — V^W an d a change of variable V((f>) = 
VqV(4>), we arrive at the action of the scalar field in the 
standard form 



'Scalar 



g^d^d^ 



exp — K4<7 



2n 

71+2 



(12) 



If the scalar field is taken to be the Higgs field H, be- 
cause the multiplicative factor exp ^—K^a^J -^p^j does 

not change the minimum of V, the Higgs vacuum expec- 
tation value, (H), is unchanged. 

For the gauge field, the action is given by 

Scaugc = J d 4 xd n y^VhG AB G CD F AC F B D, 

where Fab is the gauge-invariant field strength tensor in 
the bulk, and g„ is the 4 + n dimensional gauge coupling. 
Again we are only limited to the zeroth KK mode, which 
is the 4D effective field tensor F^. Carrying out the 
transformation Eq. [5J we have 



Sea 



1 



d x \j — g exp n^a 



2n 
n + 2 



(14) 



where we have defined F M „ = VqF^. Note that the 4D 
effective gauge couplings gl 



gl = 9l exp 



2n 



(15) 



is radion dependent. 



Similar techniques can be applied to the zero mode of 
the Dirac field ip with mass m: 

Sbirac = / d 4 xd n yV=gVh (-i^TiD^ - raw) , 

(16) 

where e v is the vielbein, I\ is the gamma matrix in 4 + n 
dimension, and is the covariant derivative. After the 
conformal transformation, we get the canonical action 



d x\f — g 



K4 

m exp I cj \ , 

11 2 V n + 2 



2n 



by redefining the field ip as 

K4 



^ = exp I — 

v \ 4 V n + 2 



2n 



(17) 



(18) 



We can read out the radion dependence of the fermion 
mass m from Eq. 1171 



m 



mexp 



K4 



2n 



(19) 



Since the fermion mass is given by y(H) in the standard 
model, and that (H) has no radion dependence, we may 

conclude that y cx exp (^—'y&^J ^fj) ■ 

In the above, we assume that electroweak symmetry is 
broken spontaneously by the condensation of the higher 
dimensional Higgs field [4l|. In [l3[, the authors con- 
struct an SM-like model from a UED model with a Higgs 
field in 5D. In particular, to generate fermions with mas- 
sive and chiral zero mode in 4D, they have to use two 
fermion fields with opposite chiralitics in addition to a 5D 
scalar field. Yet the radion dependence of the fermions in 
their model is the same as that we obtain from the above 
simple analysis. In the literature, there is another class of 
models, called Higgsless models, in which the zero modes 
of the KK tower of the fermions are removed by choosing 
suitable conditions at the extra dimensions rather than 
by the Higgs mechanism. For the Higgsless model consid- 
ered in [l3j, the radion dependence of the fermion mass 
is also the same as Eq. [T5] However, we should men- 
tion that in other types of models, e.g. the "gauge-Higgs 
model" in 42], the Higgs VEV may depend on the ra- 
dion, and thus the radion dependence of the fermion may 
be different. 



III. VARIATION OF CONSTANTS IN UED 
MODELS AND ITS IMPACTS ON THE CMB 
POWER SPECTRUM AND SNE IA 

In the last section, we have formally derived the ra- 
dion dependence of the gauge couplings and the fermion 



4 



masses. Now we are in a position to discuss the radion de- 
pendence of the constants, with the understanding that 
we are interested in the epoch of CMB formation and 
afterwards. Here we shall first discuss the radion depen- 
dences of a, m e , and Aqcd in Sect. lIII A~l We then discuss 
in Sect. lIII Bl thc effects of variation of these constants on 
the CMB power spectrum and on SNe la respectively. 



A. Variation of the relevant constants in the UED 
models 

Recall that we work in the Einstein frame, and so G 
is constant. For CMB, the relevant constants are the 
electron mass m e , electronic charge e and proton mass 
m p , whereas m p is important for SNe la. 

Among the fermion masses affected by the variable 
radion, the variation of m e is most relevant during the 
epoch of CMB formation. Since the quark mass contri- 
bution to the proton mass (or neutron mass) is only a few 
per cents or so, we neglect the effect of variation of the 
radion on nucleon masses. For notational convenience, 
we use the subscripts cmb and o to denote the values 
at CMB recombination and today respectively, and we 
define 



-1/2 
PCMB 



exp 



K4 



-OCMB 



2n 



(20) 



Note that Eq. [10] implies <tq 
we have 



and so po = 1, and thus 



-1/2 

"*e,CMB = PcMB TOe :0- 



(21) 



From Eq. 1151 we now derive the radion dependence of 
the fine structure constant a and the QCD scale parame- 
ter Aqcd- In SM, it is conventional to define on = gf /An, 
i=l, 2, 3, to denote the U(l), SU(2), and SU(3) coupling 
constants respectively, and so we have on oc p . Using 
the usual SM relation a -1 = a^ 1 + , we deduce the 
radion dependence of «cmb : 



aCMB — "OPcMB- 



(22) 



The variation of the strong coupling a 3 induces vari- 
ation of the QCD scale parameter Aqcd- As in usual 
QCD, the running of a 3 from some high energy scale Q 
to some lower energy scale \x is given by the one-loop 
renormalization group equation 



a 3 (Q 2 ) 



as(M 2 ) 



1 



■(33- 2n/)ln^ 



(23) 



where nj is the number of quark flavors with mass less 
than Q. When Q is sufficiently low, the effective coupling 
becomes very large, and this scale, Aqcd> is defined as 



A, 



QCD = M ex P 



-67T 



L(33-2n / )a 3 (M 2 ) 



(24) 



Making use of Eq. [24] and matching the coupling at the 
thresholds m c and mj, Aqcd can be expressed in term 
of az(Mz) at the predetermined energy scale of Z boson 
il 



A 



QCD 



M 5 



V Ml 



2/27 



exp 



-2vr 



9a 3 (M z ) 



(25) 



Thus we can relate Aqcd, cmb to its present value 
Aqcd : q: 



A 



QCD, CMB 



Aqcd.oPcmb 7 ex P 



2tt 



9a 3 ,o(M z ) 



(1 — Pcmb) 



Because the contribution from Aqcd dominates the nu- 
cleon mass, we simply take the nucleon mass to be pro- 
portional to Aqcd ■ Variation of m p will affect both the 
CMB power spectrum and the SNe la. 



B. Effects of variation of the constants on the 
CMB power spectrum and SNe la 

In this subsection, we study how the variation of the 
constants discussed in Sect. IIII Al affect the CMB power 
spectrum and SNe la. 

Since the effects of variation of a [20I l2ll I22L [23l , m e 
[24l [25l ] individually, and a and m e together [3l[ have 
been discussed quite extensively in the literature, we shall 
briefly discuss their impacts only. The ionization history 
has been modeled accurately by the code RECFAST, in 
which the ionization fractions of hydrogen and helium, 
and the matter temperature Tm are evolved (see 44[ and 
references therein) . The variation of a and m e modify the 
ionization history by varying the binding energies of hy- 
drogen Hi and helium Hen, Thomson cross-section cry, 
the recombination coefficients 1Z and the two-photon de- 
cay rates of hydrogen and helium. The binding energies 
scale as a 2 m e ; <jt is proportional to a 2 m~ 2 , and the two- 
photon decay rates vary as a 8 m e [HI]. In the literature, 
the recombination coefficients for hydrogen and helium 
are usually parametrized as a function of matter temper- 
ature. The recombination coefficient 1Z can be expressed 
as [13 



K = ^8vr(2/+l) 



JfeT; 



A I 



27T771, 



exp 



( B n 



M 



v n iy 2 dy 
ev-1 ' 



(27) 

where B n is the binding energy of the nth state and a n i is 
the ionization cross-section. The asterisk indicates that 
the sum should be regulated. Using the fact that B n 
scales with m e a 2 and the ionization cross-section varies 
as m7 2 ( 



[46( , we can derive [2l|, [47 



on 

da 
dTZ 



a \ oTm 



— — [ 2R 

m R 



T dU \ 
TM df^) 



(28) 
(29) 



5 



which enable us to get the radion dependence of 1Z using 
the empirical fitting formula for 1Z. 

Among these effects, the change in the binding energy 
of hydrogen, _B„, is the most significant, and so the qual- 
itive changes in the power spectrum can be understood 
based on the change in B n . 

The variation of Aqcd causes the hadron masses to 
vary. This will cause the energy density of the uni- 
verse to change, and hence the expansion rate. A con- 
sistent theory should address where the energy goes to 
or comes from. It may arise from interaction of the SM 
particles with the radion. But without being too model- 
dependent, we simply ignore the soure/sink of the baryon 
energy density due to the variation of Aq CD . A marked 
effect of the variation of the baryon energy density is the 
change in the ratio of baryon to photon energy density, 
usually defined as 



-Rcmb 



3pb,CMB 



4p 7 ,CMB 

which is related to the sound speed c s as 



<?.= 



1 



3(1 + -Rcmb) 



(30) 



(31) 



The effects of baryon on the power spectrum through 
-Rcmb are as follows [H]]. The presence of baryons re- 
duces c s and hence the sound horizon; more remark- 
ably, it shifts the equilibrium point of oscillations, which 
breaks the symmetry between the odd and even peaks. 

In Fig. [TJ we show the power spectrum with /?cmb = 1 
and Pcmb = 1-05. We also compare the effects of indi- 
vidual variations of m e , a and Aqcd for pcmb = 1-05 on 
the power spectrum in the figure. Note that for Pcmb = 
1.05, m eiC MB/«i e ,o, acMB/ao and Aq C d,cmb/Aq C d,o 
are 0.98, 0.95 and 0.74 respectively. As the change in 
m e is small, we only see a small shift to lower I scales 
due to later recombination caused by the smaller bind- 
ing energy. Because the change in a and its weight on 
the binding energy of hydrogen are greater than those of 
m e , the shifting of the peaks to small I scales is more 
pronounced compared to changing m e alone. Later re- 
combination implies that the visibility function is broad- 
ened, and so the recombination takes a longer time to 
finish and the damping effect becomes more significant 
in the large I scales. We also expect the early integrated 
Sachs- Wolfe effect to be reduced for having less residual 
radiation during a later recombination. 

Because of the exponential factor in Eq. [26j a large 
change in Aqcd results, which means that the effect on 
the CMB power spectrum is appreciable. There are two 
major effects arising from the variation of Aqcd- First, 
a decrease in baryon density reduces the Hubble rate, 
and this implies that photon decoupling takes place later. 
The consequences of later recombination are that the 
sound horizon is larger, which corresponds to shifting of 
the peaks to smaller I scales, and the damping in high 
I scales is reduced. Second, a decrease in baryon mass 



3. 




FIG. 1: The CMB temperature power spectrum with pcmb = 
1.05, and the individual contributions to the power spectrum 
due to the variations of m e , ot and Aqcd are plotted against 
the one with pcmb = 1. See the text for discussions. 



causes the sound speed to increase, which not only results 
in a larger sound horizon, but also shifts the equilibrium 
point less. The latter helps to restore the symmetry be- 
tween the odd and even peaks; in particular, the power 
of the first peak is reduced whereas the power of the sec- 
ond peak is boosted. Indeed, we see that peaks shift to 
smaller I scales due to the larger sound horizon caused 
by these two effects. 

The formation of CMB is basically determined by the 
values of e, m e and m p at the last scattering surface, and 
so the CMB constraint limits the variation of these con- 
stants during the CMB recombination epoch, z ~ 1100. 

We now turn to discuss variation of the constants that 
are relevant for SNe la. SNe la are believed to be ex- 
plosions of white dwarfs that have approached the Chan- 
drasekhar mass [49| . Mch oc G~ 3 / 2 mp 2 , which depends 
on the fundamental constants and is not sensitive to the 
structural details of the white dwarfs. Simple analyti- 
cal formula for light curves of SNe la [5(| predicts that 
the peak luminosity is proportional to the mass of nickel 
produced, which in turn is approximately proportional 
to M C h- In HI, the authors constrained the variation of 
the gravitational constant G using SNe la data, and it 
was found that the effects of of varying G on SNe la light 
curves are basically captured by the Chandrasekhar limit 
after comparing the results of simple Mch consideration 
with full numerical computations. Similarly, varying the 
nucleon mass causes Mch, and hence the intrinsic lumi- 
nosity of the SNe la to evolve. Following [34], [35|, we 
shall constrain the variation of Aqcd using the Mch ar- 
gument. Supposing that the luminosity L cx m~ 2 , we 
have 



M = M + 5 log 



(32) 



where M is the absolute magnitude and the subscript 
denotes the local value. 
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Since there is no generic theory governing the evolution 
of the radion, in this paper, we shall parametrize the 
evolution of psNe as a linear function of redshift z: 

PSNc = 1 + (33) 
Hence the apparent magnitude m can be written as 

TO = M + 51og-^- + 51ogd L + 25, (34) 
rripfi 

where c?l is the luminosity distance, and m p oc Aqcd.snc 
with Aqcd.sno evaluated as in Eg. 121)1 




IV. NUMERICAL CONSTRAINTS ON THE 
EVOLUTION OF THE RADION BY CMB DATA 
AND SNE IA AND DISCUSSIONS 

In this section, we shall constrain the radion evolution 
using the three-year WMAP data [H and then the SNe 
la data [Hlj . 

Since many free parameters are involved in CMB cal- 
culations, we make use of the Markov Chain Monte Carlo 
(MCMC) method implemented by the engine CosmoMC 
[52 ] , which searches for the region in the parameter space 
that maximizes the likelihood function, to derive the con- 
straints on the free parameters. The theoretical CMB 
spectra are calculated by the Boltzmann code CMB- 
FAST "5^ . We vary the following set of parameters: the 
Hubble parameter, iJo, the baryon density, u>i, = fl^h 2 
(h = i?o/100 km s -1 Mpc" 1 ), the cold dark matter den- 
sity, oj c = Q, c h 2 , the reionization redshift, z Te , the pri- 
mordial fluctuation amplitude, A s , the spectral index, 
n a , and pcmb- Flatness of the universe is assumed in all 
the calculations. 

Assuming a top hat prior of [64,80] km s" 1 Mpc -1 for 
H , we get the 95% confidence interval (C. I.) for pcmb 
to be [0.980, 1.021]. However, we find that Ho is not 
constrained well within the assumed interval, and, more 
importantly, Hq exhibits strong degeneracy with pcmb- 
To tighten the constraint on pcmb, we further fix Ho to 



be the HST Key Project result, 72 km s" 1 Mpc" 1 [54] 
to break the degeneracy. The 95% C. I. for pcmb now 
reduces to [0.991,1.012], and the marginal distributions of 
the parameters are shown in Fig. [21 where we see that the 
standard cosmological parameters are constrained within 
the usual regions. 

In Ref. [25j, the radion evolution in the brane world 
scenario, in which only m e is expected to be radion de- 
pendent, is constrained by CMB data. Compared to that 
in the brane world model, the constraint on pcmb in the 
UED scenario is tightened by about a factor of 2, thanks 
to the fact that a and Aqcd, in addition to m e , also vary. 

Since the upcoming Planck satellite mission is going to 
measure the temperature power spectrum to as high as 
I ~ 2500 and the i?-polarization spectrum to I ~ 1500, we 
expect that there will be tremendous improvement in the 
constraint on pcmb- We can forecast the improvement 



FIG. 2: The marginal distributions of the free parameters, 
constrained by the three-year WMAP data. In addition to 
the free parameters ujb = Qi,h 2 , uj c = Q c h 2 , z rc , n B , A B and 
Pcmb , shown also are the derived distributions of the density 
parameters of matter (f2 m ) and cosmological constant (Qa). 
The Hubble parameter is fixed to be 72 km s" 1 Mpc -1 . The 
maxima of the distributions are arbitrarily normalized to 1. 



that Planck will bring quantitatively using the Fisher ma- 
trix, which has been widely used to predict the expected 
uncertainties in future experiments (see e.g. (55j). Under 
the assumption of Gaussian perturbations and Gaussian 
noise, the Fisher matrix takes the form 



^-s ^ dCx ; ,„ s_i<9Cy/ 



X.Y 



(35) 



where pi is the ith free parameter and Cxi is the Ith 
multipole of the observed spectrum of type X, which 
can be the temperature, temperature-polarization and 
i?-polarization spectra. The experimental precision is 
encoded in the covariant matrix Covixy- Using the ex- 
pected precision of the Planck satellite, we find that the 
current bound on pcmb will be tightened by a factor of 
8 or so when Hq is fixed or free. 

Now we turn to the constraint on the radion evolution 
by SNe la data. To be consistent with CMB calculations, 
we also assume flatness of the spatial curvature in evalu- 
ating d,L in Eq. 1341 Ho is analytically marginalized over 
[561 ] . Assuming a linear variation of psno as in Eq. [331 we 
get the joint x 2 distribution of the density parameter of 
matter, fi m , and £, as shown in Fig. [3] Correspondingly, 
the 95% C. I. for C is [-0.025, 0.051]. 

To tighten the constraint on f , we impose a Gaussian 
prior on fl m using the WMAP plus SDSS constraint on 

n. 



2fifi+ 026 



57] 



We simply take the mean of the 
lower and upper error bars as 1 a for O m . The resultant 
95% C. I. for after numerically marginalizing over f2 m , 
is [-0.005, 0.015], or [0.998, 1.006] for p SNc . The bounds 
on ( are asymmetric due to the fact that the value of Q m 
derived from the SNe la gold sample [5l| is higher than 
that favoured by WMAP and SDSS data (see [57|). 
Because the constraint on psNe is solely due to the 
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FIG. 3: The joint \ 2 distribution of Q m and f , constrained by 
SNe la data. The dashed and solid lines are the 68% and 95% 
condidence contours respectively. The horizontal (black) lines 
enclose the 1 a region of the Gaussian prior that we impose 
on Q m , while the vertical (red) lines entail the 95 % C. I. on 
£ after marginalizing over Q. m . 



[H, [Hill or A QCD . 

More drastically, variation of constants may even mim- 
ick the effect of dark energy, but this is increasingly un- 
tenable as more and more evidences, such as the CMB, 
necessitate the existence of the cosmological constant or 
suchlikes. 

Our argument of the effect of variation of Aqcd on the 
light curves of SNe la rests on the fundamental Chan- 
drasekhar mass; nonetheless it is just the first step and 
should be checked by detailed simulations that the fine 
details of the explosion per se do not give rise to signifi- 
cant corrections. This may be a daunting task given that 
numerous nuclear reaction cross-sections can be affected. 
See [gH for a discussion on the effect of varying strong 
coupling constant on the fusion reaction in the SNe la. 

Finally, we would like to compare the constraints on 
p z at different redshifts z. The constraints from BBN 
[3, CMB and SNe la are displayed in Table Q] It is 
meaningful to consider p/p = (1 — p z )/\At\. For SNe 
la, we use the mean redshift of the gold sample, 0.4, 
to calculate At. For the constraints on Aqcd (Eq. ISM)) , 
they can be translated to the bounds on Aqcd/ Aqcd, o : 
[—1.4,2.8] x 1CT 11 yr _1 . Again, we mention that the 
simplest light radion models have been ruled out by test 
of general relativity and fifth force experiments although 
we don't include them here. 



variation of Aqcd , it is useful to translate the constraint 
on psNe to the constraint on Aqcd: 

Aqcd.sno = Aq C d,o(1 + (36) 

with the 95% C. I. for £ to be [-0.30, 0.15]. We stress 
that Eq. alone can be interpreted as a constraint on 
the variation of Aqcd , independent of UED model. 

It is worth mentioning that we have also considered 
parametrizing the variation of Aqcd as Eq. [36] and con- 
fronting it with the SNe la data directly, and the resul- 
tant constraint on £ is similar to that derived from the 
constraint on £. In other words, the bounds on Aqcd are 
not very sensitive to the two type of parametrizations, 
Eq. [32] and Eq. [Ml 

Workers in SNe la community have been aware of the 
possibility of evolution of SNe la light curves. The light 
curves are usually adjusted based on the Phillips relation 
[Bil l or other alternative schemes such as the stretch fac- 
tor [5f|. These schemes implicitly assume that there is 
one underlying parameter governing the light curves, al- 
though the nature of this parameter has not been clearly 
identified. However, scatter about the mean relation al- 
ways exists, and a recent study exacerbates the inade- 
quacy of the one-parameter description [60]. Although 
most theorists suspect that another underlying param- 
eter is related to the metallicity Z of the progenitor 
[6lT l62l . HH, observations so far do not reveal any defi- 
nite Z dependence [H, [65[ . So far most of the studies 
concentrate on metallicities or the two-component model 
[H, [13] , we would like to emphasize that the second pa- 
rameter may be due to variation of constants, such as G 



TABLE I: The constraints on p from BBN, CMB and SNe 

la. 

Redshift z I Observations 1 95 % C. I. for p z |95 % C. I. for p/p 

(i<r 13 yr- 1 ) 

^ 10 10 BBN [18] [0.9992, 1.0027] [-2.0, 0.59] 

~ 1100 CMB [0.991, 1.012] [-8.8, 6.6] 

~ 1 SNe la [0.998, 1.006] [-14 , 4.7] 



V. CONCLUSION 

In this paper, we consider variation of fundamental 
constants in the context of Universal Extra Dimension 
models, where not only the gravity, but also the stan- 
dard model fields can propagate in the extra dimensions. 
We explicitly show the radion dependence of various con- 
stants by going through the dimensional reduction pro- 
cedure. It turns out that the constants relevant for CMB 
power spectrum are a, m e and Aqcdj while it is Aqcd 
that is important for SNe la light curves. The effects of 
the variation of the constants on the CMB power spec- 
trum and SNe la are discussed. For CMB, we have 
done detailed numerical calculations to get the power 
spectrum. The effect of variation of Aqcd on SNe la 
light curves is estimated based on the fundamental Chan- 
drasekhar mass. We then use the three- year WMAP data 
to constrain the radion evolution at z ~ 1100, and the 2 er 
constraint on p/p is [-8.8, 6.6] x 10" 13 yr" 1 . The SNe la 
gold sample yields a constraint on p/ po, at redshift from 
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to 1, to be [-14,4.7] x 10~ 13 yr^ 1 . Furthermore, the 
constraints from SNe la can be interpreted as bounds on 
Aqcd/Aq C d,o, [-1.4,2.8] x 1CT 11 yr" 1 . 
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